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Abstract. We study non-degenerate (Petrov type I) silent universes in the presence 
of a non-vanishing cosmological constant A. In contrast to the A = case, for which 
the orthogonally spatially homogeneous Bianchi type I metrics most likely are the only 
admissible metrics, solutions are shown to exist when A > 0. The general solution is 
presented for the case where one of the eigenvalues of the expansion tensor is 0. 

PACS numbers: 0420,0440 
1. Introduction 

The electric and magnetic part of the Weyl tensor with respect to some unit timehke 
vector field u°' are defined by 

with rjabcd the 4-dimensional Levi-Civita tensor. Silent universes are defined ^2] as 
irrotational expanding dust models, which are solutions of the Einstein field equations 

Rab - ^Rgab + J^gab = Tab = fJ'UaUb (2) 

(wherein /i 7^ 0) with the additional property that the Weyl tensor is purely electric 
with respect to the fiuid 4- velocity 

Hab = 0. (3) 

The Petrov type is then necessarily O, D or I [T^, with the type O solutions given by 
the Friedman-Robertson- Walker dust metrics. All Petrov type D solutions are known 
explicitly too: they are characterised P by the fact that the Weyl tensor is degenerate 
in the same plane as the shear and are given by the Szekeres ^H] family, including the 
subcase of the Ellis jUj LRS class II dust models and such well known examples as the 
Lemaitre-Tolman-Bondi model and the orthogonally spatially homogeneous Kantowski- 
Sachs model. These models can be easily generalized to take into account a non-zero 
cosmological constant [51 IT^ 113]. 
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For the general Petrov type I case the situation is rather different, the only known 
silent universe being the orthogonally spatially homogeneous Bianchi I model. The issue 
of the integrability conditions of the inhomogeneous silent universes of Petrov type I has 
been a source of confusion and discussion during the last decade El UHl HH HI] and, 
although the field equations of the silent models are consistent when linearised about 
a Friedman- Robert son- Walker background, the full non-linear consistency problem re- 
mains — to this day — unsolved. A detailed study of the constraints so far only has 
led to the conjecture that the Bianchi type I models are the unique Petrov type I silent 
universes. 

At first sight it looks as if the presence of a cosmological constant is of little 
relevance in the consistency problem (it might well be that this is indeed so in the 
generic situation). However, one of the special cases where one can show that the 
constraints with A = are inconsistent (namely the case where one of the eigenvalues 
of the expansion tensor is 0), surprisingly turns out to admit explicit type I solutions 
when A > 0. 

2. Basic equations 

We first interpret the cosmological constant as a constant negative pressure p = —A and 
rewrite the field equations (jSJ in the form (also substituting fi hj + p ^ 0) 

Rab - ^Rgab = (/i + p)UaUb + PQab (4) 

and we further follow the notations and conventions of the Ellis-MacCallum orthonormal 
tetrad formalism 9j. Herein, the normalized 4- velocity u"" plays the role of the timelike 
basis vectorfield Cq of the tetrad {cq, e'^}. Greek indices take the values 1, 2 and 3 and 
refer to tetrad components with respect to e^. Consequently, for spatial tensorfields (i.e. 
for which all contractions with u"' vanish) only the Greek indexed components are not a 
priori zero. For spatial vectorfields Va and rank 2 tensorfields Sab these components Va 
and Sai3 form (3 x 1) resp. (3 x 3) matrices denoted by boldface symbols v and S and we 
will write Sa for the eigenvalues of S. As a trivial example, the tensor hat = gab + UaUb, 
which projects orthogonally to u"", satisfies h = 1, i.e. ha/3 = ^a/s- 

The basic variables in this formalism are 24 independent linear combinations of the 
Ricci rotation coefficients or, equivalently, of the commutator coefficients 7^^ associated 
with the tetrad; these are the objects (a < P, with additionally n/3a = Uap) and 
defined by 7^^ = e^sji''^ + S'^a^ — d^as, the components VL^ of the angular velocity of the 
triad Ga with respect to the local 'inertial compass' and the components Mq,, uja and O^p 
of the spatial kinematic quantities (acceleration vector, vorticity vector and expansion 
tensor respectively). We also use the decomposition = + 1/365q,/3, with aab the 
trace-free shear tensor and the expansion scalar. It turns out that there is a slight 
notational advantage in using instead of the variables n^p (a 7^ /3) and the variables 
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qa and defined by 

Here and below these expressions have to be read modulo 3, so for example a = 3 gives 
ni2 = 2^3 + rg. 

The basic equations in the formalism are the Jacobi-equations (77)-(81) and the field 
equations (82)- (84) of jH] (where the print error in (83) was corrected in [d) together 
with the energy-momentum conservation equations or 'contracted' Bianchi identities. 
Below, the labels (77)-(84) refer to the corresponding equations in [H]. 



For problems in which conditions on E and H are imposed one can proceed as 
follows. Both tensorfields being spatial, the only not a priori vanishing components are 

Ea(3 = CaOpO and HaiS = -Vao"^^Cmnl30 (5) 

and can be expressed in terms of the basic variables and their derivatives in two ways. On 
the one hand, one can take the definition of the Weyl tensor through the decomposition 
of the curvature tensor (see (3.45) or (3.50) of T|), using the formula for the curvature 
components in terms of the Ricci rotation coefficients and their derivatives (see the 
appendix) and using the expressions Rafs for the Ricci tensor components originating 
from the field equations (jH), to yield expressions which we denote by E^^, H^^, on the 
other hand one can deduce [0] expressions -E^^, H'^^ from the Ricci identities for m", 
which are given in 1 + 3 covariant form by (A3) and (A6) of ^Uj. Both ways are of 
course equivalent via the Einstein and Jacobi equations: one can check that 
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Ka - H^a = ^ ^[(81)"" - (81)°+^"+^ - (81)°-^"-^] + ^(77), 

K+ia-i - = ^ ^(80)" + (81)"+i"-i - ^(83)" 

- H^-ia+i ^ -^(80)" + (81)"+i""i + i(83)^ 

where for example (80)^ and (81)^^ j^g^n '(80) with 7 = 3' and '(81) with a= l,p = 2' 
respectively. 

These components are treated as extra variables and form matrices E and H, the 
symmetry and tracelessness of which is automatically fulfilled when using E'^ and H'", 
but for E'-^ and H"-^ only follows from specific Einstein and Jacobi equations — as is 
evident from ©. 

Finally, the 16 remaining Bianchi identities are added to the basic equations; 
they can most easily be translated in terms of the above variables from their 1+3 
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covariant forms (e.g. (A7)-(A12) of [TU]) by means of the connecting formulas presented 
in the appendix. Below, the labels (Ax) refer to the corresponding equations in pUj . 

Returning to the case of silent universes, we have to cope with the basic equations 
under the additional conditions 

H = 0, a; = 0, u = 0, e^O, (^i-E2)(^2-^3)(^3-^i) 7^0,(7) 

the last condition expressing the fact that we restrict ourselves to the Petrov type I case; 
the third condition follows directly from the energy-momentum conservation equation 
(A8), while (A7) determines the evolution of the matter density /i along the matter flow 
lines, 

fM = doi2 = -{n + p)&. (8) 

We next recall some basic results from J]. Under the conditions (jZj) the divH 
Bianchi identities (A12) become algebraic, 

eafsy'Es'' = 

{eafB'y denotes the usual 3-dimensional Levi-Civita skew symbol) and expresses exactly 
that the antisymmetric part of a • E vanishes, i.e. that a and E commute. Therefore 
it becomes advantageous to work in an orthonormal eigenframe of E, which thus is 
simultaneously an eigenframe of a and 6: this fixes the triad {e^} completely since the 
Petrov type is I and drastically decreases the number of variables. First, the sets {Eais} 
and {Oa/s} are reduced to the respective eigenvalue sets {Ea} and {9a}- Second, the 
off-diagonal part of (A9) 

Eap = doEap + 2n'''eyS{aEp)^ = -QEap + 3a'^{aEp)y -^if^ + P)^«/3 (9) 

shows that an eigenframe of E and a is also Fermi-propagated along the matter fiowlines, 
i.e. = 0. Third, as shown in [T] for Petrov type I fields, the diagonal parts of (A6) 
and the H Bianchi identities (AlO), 

= H^^ = (curla)«„, = if^ = -(curlE)„„ (10) 

(where curl is the covariant spatial curl operator, the definition of which is also given 
in PO]), are equivalent with the fact that the Weyl principal vectors are hypersurface 
orthogonal, i.e with the vanishing of all Ricci rotation coefficients j^-i a+i 
formalism, with n^a = 0. More precisely, Barnes shows that ()10|) implies either this fact 
or else that the fluid is shear free (a = 0) and thus [2] static, the latter possibility being 
contradictory to the condition O 7^ 0). 

Under the previous conditions the Jacobi equations (77), (79) and the diagonal part 
of (81) are identically fulfllled, so that the the diagonal components of and H*-^ are 
the same (see ©) and vanish identically. There are neither extra calculations involved 
in expressing H^^^ = 0, which implies the evolution equations 

dor a = -Oa+lTa ^^^^ 
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The expression for Ei together with the (0 0) field equation (82) give 

do9^ = -E^-9l-^{fi + 3p) ^^^^ 
E33 = —{Ell + -^22)- 

Substituting the expressions (fTT|) in the Jacobi equations (80)° and (81)"+^"-\ these 
readily determine the spatial derivatives daOa+i and daOa-i to be 

daOa+l = '2ra{9a+l — 9a) ^^^^ 

Under (fTT|) and (fTSjl the (0 a) field equations (83) become identities, see (0). The 
(a + 1 a — 1) field equations together with the only remaining Jacobi equations (78)" 
imply 



(14) 



daQa-l = 2ra{qa-l + ^q-i) 

while, with (fT^ . the (aa) field equations yield 

-6a0a+i — ij^^a-i + -/ 

At this point all Jacobi and Einstein equations have been dealt with. Finally, the 
diagonal part of the E Bianchi identities (A9) yields 

d^Ea = {da - e)Ea + Ea+lOa-l + ^a-l^.+l + + p){0 " 3^J, (16) 



while the divE Bianchi identities (All), and then the off-diagonal part of (AlO), lead 
to the following expressions, where the variables ma represent the spatial gradients of 
the matter density, nia = dafi- 

daEa = -^fna — 2ra{Ea+l — Ea) — 2qa{Ea — Ea-l) 

daEa+1 = -\ma + 2r„(Eo+i - Ea) (17) 
o 

daEa-1 = —-^^a — 2qa{Ea~l — Ea) ■ 

o 

Herewith all basic equations are satisfied identically. 



3. Integrability conditions 

Propagating the expressions for the spatial gradients of Ea along the matter fiow lines, 
and using the commutator relations (75) of [Oj and the above deduced equations, one 
obtains 

{Ea-l — Ea+l)daOa = 6Ea[qa{0a-l " ^a) + ^q(^q+1 " ^a)] + -^^ai^^a+l " ^a-l) (18) 
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and 

{Ea-i - Ea+i)doma = 4(/i + p)[qo,{Eo, - Ea^i){6a - 9a~l) + ra{Ea - -Ea+l)(6'„ - 6*0+1)] 
+im„(/i + p)(0o_i -m„(0 + 0„)(Eo_i -E^+i) (19) 

Further integrability conditions can be obtained by considering the [da-i, da+i\Ea com- 
mutator relations, but this is at present not our concern. 

As Barnes noticed in jT], the expansion tensor of type I silent universes is non- 
degenerate: the time evolution of, for example, 9i = 62 would immediately imply 
El = E2, see (fT^ . Looking at the equations (jSjl. lfT^ and (fTBj) one sees that do is a deriva- 
tion on the polynomial algebra M.[Ei, E2, 61, 02,03, fj]- Therefore it is fairly straightfor- 
ward to investigate whether a certain algebraic relation between the 6 mentioned vari- 
ables is consistent: to obtain new algebraic relations one can take time derivatives of a 
newly generated algebraic relation or calculate the resultant of it and an already deduced 
algebraic relation (with respect to an appropriate variable which thereby is eliminated 
from these two). From the above, the most natural Ansatze seem to be ^1 = or Ei = 0. 
However, for the last one we get the same inconsistency result as for general vacua (cf. 

Proposition There exist no silent universes (with or without cosmological con- 
stant) for which one of the eigenvalues of the Weyl tensor vanishes. 

Sketch of proof: We denote the resultant of two polynomial ex- 
pressions Pi[xi, . . . ,Xn] and P2[xi, . . . ,Xn] with respect to the variable Xi by 
reSa;. (Pi, P2)[xi, . . . ,Xi, . . . ,Xn] (where Xi means that the variable Xi is supprimed). Sup- 
pose e.g. El = 0. From (fTHj) with a = 1 we get 

Pi[E2, 01, 02, 03, /i] := ^(-2^1 + 02 + 03) + (03 - 02)E2 = 0. (20) 

Put P2[0i,02,03,t^] := 18TeSE,{doPl,Pl)[0l,02,03,lA/il^ + P) and P^[0i,02,03, l^] ■■ = 

TesE2{doP2, Pi)[0i, 02, 03, tA/ {02-03)- Herewith the resultant res^(P2, P^)[0i, 02, 03\ yields 

3(^2 - 03)\'^0l -02 - 03) {<O0l - 3^2 - 3^3 - 7^102 " T^i^s + 14^2^3) = 0.(21) 

The first possibility to satisfy this relation, 03 = 02, is not allowed because of the non- 
degeneracy of the expansion tensor; the second possibility 20i — 02 — 03 = comes down 
to the first one, since expressing 9o(26'i — 02 — 03) and substituting 0i = {02 + 6*3) /2 
yields —{02 — ^3)^/2, from which we get 

^"4(^1, 02, 03] ■■= G0l - 30l - 30l - 70102 - 70103 + 14^2^3 = 0. (22) 

Here, P4 is homogeneous and quadratic in 6*1, 6*2 and 03, but not positive definite. 
Therefore we proceed and calculate successively 

TeSE,{doPA,Pl)[0l,02,03,fi] ■■= {02 - 03)Ql[01,02,03,f^]/Q 
TeS^{Ql,P2)[01,02,03] ■={20l- 02 - 03)Q2[01,02,03] 
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and get 

res,3(Q2,P4)[^i,^2,^3) := -96(^i - ^a)'^!^!, ^2], (24) 

where P5 is a polynomial of total degree 4 in 9i and 02- Proceeding in exactly the same 
way for P5 as for P4, one gets polynomials Pi [6*1, 6*2, 6*3, /i] and P2 [6*1, 6^2, 6^3] . Finally, 
elimination of ^3 and 62 from P2, making use of P4 and P5 respectively, results in a 
polynomial of degree 40 in ^1, with coefficients depending on p = —A. This implies 61 
is constant and hence, by P5, also 62 is constant. But then is constant by (fT^ (for 
a = 1) and thus B = from (jH}, which is the desired contradiction. 

Most likely, extended classes of the previous Ansatz (e.g. Ei = aE2, a constant) 
can be proven to be inconsistent by the same technique, but we do not consider these 
here. 



For the other case ^1 = we can start as above, but now we get a qualitatively 
different picture. 

From doOi = one obtains 

Pi = -^(/i + 3p), (25) 


two further time derivatives of which yield consecutively 
and 

{fi + 3p)^ + Ap{e2- 6^)^ = 0. (27) 

For p ^ (i.e. the case of vanishing or negative cosmological constant), this leads to 
fi + 3p = and thus to the inconsistency Ei = E2 = 0. Taking a further time derivative 
of (|27|). shows that it is conserved along the matter flow lines, suggesting to have a closer 
look at the possible existence of solutions for p < 0, i.e. A > 0. 



4. Solutions with = and A > 

With ^1=0 one deduces from (fT^ that ^2 = ^3 = 0. Herewith di9i = gives an 
algebraic expression for mi, which turns out to be conserved along the matter flowlines 
under the constraint (P7j) : 

mi = 2(/i + 3p) / \ (28) 

O2 — ^3 

Whereas no new information enters by propagating the constraint (jTfj) along the di 
integral curves, propagation along 82 and ^3 shows that also Qz = ^2 = 0. Herewith the 
expressions (ITH|l for a = 2 and a = 3 result in 
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and 



11 + p- 2626^ - Sqm = 0. 

It is advantageous now to introduce a constant A 
allowing one to rewrite (P7j) as 



Herewith flHUj) becomes 

4giri + (^2 + A)(03-A) =0 
and the remaining equations for the variables qi,ri,62, 0^ and /i simplify to 



(30) 

\/A), defined by p = — and 

(31) 

(32) 



doqi 

diQi 
d2qi 



-2g2 - 1a(A - ^3 







don = -02^ 

d,n = 2t\ + ^A(A + 02) 



(33) 





d^rx = 

-203^1 

0. 





^3) 



(34) 



It is clear that gi = iff 6'3 = A, while ri = iff 6*2 = — A. It is easy to check that 
the integrability conditions for the remaining variables are then identically satisfied. 
When qiTi 7^ the integrability conditions for the equations IMj) imply further 
relations, which allow one to close the system by introducing new variables u = 6'^^c?2^2 
and V = O^'^d^Os, together with their higher order derivatives u92u' = d2U, v9^v' = d^v 



and u''^62u" = d2u', v'''6sv" = d^v'. The following extra equations can then be generated: 



y2 



^391 



{02 + \)V 



4ri 



01 



2 Am 



-2Am- 



ri 



A + ^2 








Am' 




< 










d2ri = 

dov 
< div 

d2V 

-Xv' 







(^3 - \)U 

4gi 
-2\v 
2Xv 





el 



X-9. 



(35) 



(36) 



(37) 



and 



dou" = diu" = dsu" = dov" = dW = d^v" = 0. (38) 

The complete set of integrability conditions for these equations implies — and is 
identically satisfied under — the relations 

d2u" = e2u'^u"), dsv" = e3v'^iv") (39) 
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Notice that the eigenvalues of E are given by 

El = -^{92 - 03), ^2 = ^(^3 + 2^2), ^3 = -^(^^2 + 2^3), (40) 

clearly showing that the Petrov type is I (unless 6*2 or 6*3 = 0; ^^2 + ^'3 = is excluded 
as this would imply = 0). 

One can easily show (using the fact that the frame is fixed and all variables are 
invariants) that in general, with rigi 7^ 7^ uv, the corresponding solutions will only 
admit a group of isometries when $ = \1' = i. In that case it follows that 

"'--^^^"^■"'-^^"^ 

and the — in general unique — Killing vector is given by 



K = -90 + ^^d2 + ^^d^. (42) 



02 + A 03 -A 
uOl vOl 

When qiTi = or = an extra Killing vector may arise (cf. below). 



Using the fact that the tetrad basisvectors are hypersurface-orthogonal [T], it is 
straightforward to construct the metric§: 

when qiTi 7^ one finds, with fi{y) and gi{z) 7^ 0, 
ds'^ = -dt^ + dx^ + (e~^* + fi{y) cos Xxfdy^ + (e^* + gi{z) sin Xxfdz'^, (43) 
with matter density given by 

^2 3/1 {y)gi (z) sin Ax cos Ax + e'^*/i (z/) cos Ax + e^^^gi {z) sin Ax — 1 ^^^^ 
(^g-Ai _|_ j-^ (,Qg (e'^* + gi{z) sin Ax) 

whereas for gi = or ri = one obtains 

ds'^ = -dt^ + dx^ + [e-^* + h{y) cos(/2(y) + \x)fdy'^ + e^^'dz^. (45) 

with matter density given by 

JM cosjMy) + Xx) + e~ ^' 

fi{y) cos{f2{y) + Ax) - e- 

The solution (j43|) admits a (unique) Killing vector when either fi or gi is a constant 

or when fi{y) ~ y^^ and gi{z) ~ z~'^. It admits a two-dimensional group of isometries 

when both fi and gi are constants. The solution ()45p obviously admits a Killing vector, 

aligned with 62 or 63 and will admit a second Killing vector when /i and /2 are constants. 

§ A coordinate transformation is used to reduce the coefRcients of e^^* to 1; when these coefScients 
are one obtains the Petrov type D subcase with 9i ^ 62 = or 9i = 63 — 0. 
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5. Conclusion 

We have shown that, in contrast to the case of zero cosmological constant, where Petrov 
type I silent dust universes appear to be restricted to the orthogonally spatially homo- 
geneous Bianchi type I , for positive cosmological constant 'large' classes of solutions 
exist (in the sense that they contain arbitrary functions of the coordinates). The key 
equation in the construction is equation ()27|). which relates the matter density to the 
eigenvalues of the expansion tensor and follows from imposing the Ansatz that one of 
the eigenvalues of the latter is 0. It clearly shows also why such an Ansatz does not 
lead to solutions when A = 0. 

The solutions could be explicitly obtained and are given by the solutions (jlSj) and 
(I45|l ■ It is an open problem whether these solutions are the only A 7^ ones of Petrov 
type I. 

6. Appendix 

In order to avoid confusion with different conventions for the Ricci rotation coefficients, 
we prefer to define the kinematical quantities via the connection 1-forms T"-b- For the 
tetrad {eg} and dual basis {e"} the latter are defined by 

T\ = V\,e'' = (e^eb^',e,*)e^ (47) 

With Tab = Qac^^b (such that T^a = —Tab) and with -u* = Cq, the components of these 
1-forms with respect to {e*^} are then given in terms of the variables appearing in [U] by 

(48) 

From this, the expression 

Rabcd = QJ^abd — abc + ^ aeJ^'' bd — ^aed^^bc + ^ abe^^ cd — dc) (49) 

and the decomposition (3.50) in 0, one can derive expressions for the components of the 
Weyl tensor (and thus of its electric and magnetic part) in these variables, as mentioned 
in the text. 

For general purpose, we give also the following connecting formulas (the covariant 
spatial derivative divergence operator div and curl operator curl appearing herein 
are defined in e.g. ^01 and are acting on a spatial vector f " and a spatial and symmetric 
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rank 2 tensor Sab)'- 
div V — d^v'^ — la^v^ 

Va = doVa + taf^^Q'^V^ 

div Sa = dpScf — Sttf^Sa^ — eafi-yU^^ Ss^ + tiS 

Sap = doSafj + 2e^5(Qr2'^5'^)'^. 
References 

[1] Barnes A and Rowlingson R R 1989 Class. Quantum Grav. 6, 949 
[2] Barnes A 1973 Gen. Rel. Grav. 4, 105 

[3] Barrow J D and Stcin-Schabcs J 1984 Phys. Lett 103A, 315 
[4] Brans C H 1975 J. Math. Phys. 16, 1008 

[5] Bruni M, Matarrese S and Pantano O 1995 Phys. Rev. Lett 74, 1916 
[6] Ellis G F R 1967 J. Math. Phys. 8, 1171 

[7] Kramer D, Stephani H, MacCalluni MAH, Hoenselaers C, Herlt E, Exact solutions of Einstein's 

field equations (Cambridge University Press, 2003) 
[8] Lesame W M, Dunsby P K S, Ellis G F R 1995 Phys. Rev. D 52, 3406 

[9] MacCallum MAH 1971 Cosmological Models from a Geometric Point of View (Cargese) Vol 6 
(New York: Gordon and Breach) p 61 
[10] Maartens R, Ellis G F R and Siklos S 1997 Class. Quantum Grav. 14, 1927 
[11] Mars M 1999 Class. Quantum Grav. 16, 3245 
[12] Matarrese S, Pantano O and Saez D 1994 Phys. Rev. D 72, 320 
[13] Sopuerta C F 1997 Phys. Rev. D 55, 5936 
[14] Szafron D A 1977 J. Math. Phys. 18, 1673 
[15] Szafron D A and Wainwright J 1977 J. Math. Phys. 18, 1688 
[16] Szekcres P 1975 Comm. Math. Phys. 41, 56 
[17] Van den Bergh N 1988 Class. Quantum Grav. 5, L169 

[18] van Elst H, Uggla C, Lesame W M and Ellis G F R 1997 Class. Quantum Grav. 14, 1151 



